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Abstract. 

We discuss the parametric instabilities of large-amplitude parallel propagating 
(— I ' Alfven waves using the 2-D PIC simulation code. First, we confirmed the results in the 

pi !■ past study [Sakai et ai, 2005] that the electrons are heated due to the modified two 

f-^ [ stream instability and that the ions are heated by the parallel propagating ion acoustic 

\-^ ' waves. However, although the past study argued that such parallel propagating 

C/3 , longitudinal waves are excited by transverse modulation of parent Alfven wave, we 

consider these waves are more likely to be generated by the usual, parallel decay 
instability. Further, we performed other simulation runs with different polarization 
of the parent Alfven waves or the different ion thermal velocity. Numerical results 
suggest that the electron heating by the modified two stream instability due to the 

^^ . large amplitude Alfven waves is unimportant with most parameter sets. 
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1. Introduction: file preparation and submission 

Large amplitude, low-frequency Alfven waves constitute one of the most essential 
elements of magnetic fluctuations in the solar corona and solar windJH |2]. Due to 
small collisionless dissipation rates, the waves can propagate long distance and efficiently 
convey such macroscopic quantities as momentum, energy, and helicity. Since loading 
of such quantities is completed when the waves damp away, it is important to examine 
how the waves can dissipate. Among various possible dissipation processes of the Alfven 
waves, parametric instabilities have been believed to be important [3], H]. Further, the 
parametric instabilities of these Alfven waves can cause the developed turbulence, which 
are believed to produce in situ the localized structures (shocks and discontinuities) in 
the solar wind[5l [6]. 

However, most past studies using kinetic simulation code mainly discussed the 
parametric instabilities among parallel propagating waves in one-dimensional system, 
although the real solar wind is three-dimensional. Actually, arguments including 
obliquely propagating waves have been developed by several authors using the fluid 
systems[7l [HI [9l [TOl [HI [12l [13]. These studies clarifled that the growth rates involving 
the oblique waves are usually much smaller than those among parallel propagating waves. 

On the other hand, Sakai et al. |ll] recently performed two-dimensional particle- in- 
cell (2-D PIC) simulations. They reported that the left-hand (LH-) polarized Alfven 
waves are dissipated through the transverse modulational instability when ion and 
electron beta is very low (= 0.01). This is in contrast to the past studies using the fluid 
systems which indicated that parametric instabilities including obhquely propagating 
daughter waves (transverse modulation or fllamentation instability) can be dominant 
only in the high (3f (the squared normalized sound velocity) plasma. One of the purposes 
of the present study is to clarify the reason of this difference on the type of parametric 
instabilities between the fluid and the PIC simulations. 

Sakai et al. [H] also reported the electron heating through the modifled two stream 
instability (MTSI) [15] which is caused by the difference in ion and electron bulk velocities 
(cross-fleld currents). In a presence of flnite amplitude Alfven waves, the cross-fleld 
currents can initially be given by the Walen relation. However, Sakai et a/. [13] only 
discussed LH- mode in their paper. Another purpose of this study is to conflrm the 
importance of the MTSI driven by Alfven waves with both right-hand (RH-) and left- 
hand (LH-) polarizations, and also by a mixture of Alfven waves (in which the Walen 
condition is not satisfled). 

The paper is organized as follows. In section 2, we present our simulation model. In 
section 3, we present our simulation results. A summary and conclusions of this paper 
are given in section 4. 



Parametric instabilities of Alfven waves: 2-D PIC simulation 

Table 1. Parameters used in simulation runs. 
Run Polarization uJo/Qi k^CA/^i Vth,i/c Ujq 
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2. Simulation setup 

We have carried out the 2-D PIC simulation |15j with the system size L^ = 1024 — 2048 
and Ly = 512 — 1024 with 100 super-particles (both ions and electrons) per cell. 
We use periodic boundary conditions in both x and y directions, with the uniform 
ambient magnetic field {bxo = 1) in the x-direction, and the homogeneous plasma (the 
zeroth order longitudinal bulk velocity u^jq = (j = i, e) and the number density 
nio = ne = const). The ion to electron mass ratio is nii/rrie = 16, and the ratio of the 
electron plasma to the cyclotron frequency is cUe/l^el = 1 (cUj and Qj are plasma and 
the cyclotron angular frequencies of electrons {j = e) and ions {j = i), which include 
the ). Alfven velocity Ca = 0.25c, thus, when the beta ratio is 0.01 for both ions and 
electrons, electron thermal velocity is Vth,ei= y 2Te/?7ie) = 0.1c, and ion thermal velocity 



is Vth,ii= j2Ti/m.i) = 0.025c. The Debye length Vth,e/^e = 1 (grid size) and the electron 
skin depth c/uJe = 10. As initial conditions, finite amplitude, monochromatic Alfven 
waves are given 

bp = boexp{i{uJot- kox)), (1) 

Ujp = Ujo exp {i {uot - kox)), (2) 

where b = by + ib^ (complex transverse magnetic field) and Uj = Ujy + iujz (complex 
transverse bulk velocity of each species), bo = 0.5, uq and ko satisfy the dispersion 
relation of the parallel propagating electromagnetic waves in the two fiuid (electron-ion) 
system [16], 



e 



(3) 



u;^ U!{UJ + Qe) uj{u + Qi)' 

and Ujo is given by the Walen relation [T7] 

fco (^0 + ^0 
We adopt the notation that the positive (negative) ojq corresponds to the RH- (LH-) 
polarized waves. In our PIC simulation the charge neutrality cannot be assumed and 
displacement current can not be neglected due to the values assigned to Ca and uje/\^e\- 
While these parameters used in the present study are different from those in the past 
studies using the hybrid simulation code|ll [HI [19], [20], our results can basically be 
explained within the framework of the Hall-MHD theory as we will see later. 
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Figure 1. Time liistories of magnetic field energy (thick; line), ion (dashed line) and 
electron (thin line) total energies in Rmi 1. 



We performed four simulation runs with different parameters as tabulated in 
Table [H The system size is L^ = 1024 and Ly = 512 in Run 1, 3 and 4, and L^ = 2048 
and Ly = 512 in Run2. We have confirmed that the results remain essentially the same 
when we extended the transverse system size to Ly = 1024. Run 1 is almost the same 
as the run in Sakai et a/.[Tl](LH- parent wave satisfying the Walen relation). In Run 2 
the RH- parent waves are given, and in Run 3, magnetic and the velocity perturbations 
are given independently (and thus the Walen relation is not satisfied) in order to model 
the fluctuations at the chromosphere. In Run 4, we introduce higher ion temperature in 
order to examine the roles of kinetic modulational instability [H [191 120]. Initial Alfven 
waves are all forward propagating in Runs 1, 2, and 4, while forward and backward 
propagating waves are mixed in Run 3. 

To discuss the propagating direction of Alfven waves, the two-dimensional (x — y) 
spacial Fourier transforming analysis were performed in this study. Since we focus on 
the low-frequency waves in this paper, the vertical and horizontal axes {ky and k^^) of 
the figures indicating the Fourier analysis (Fig. [3], [5l [Tj) are normalized by the ion skin 
depth {c/ui = 40). 



3. Numerical results 

First, we show the results of Run 1, in which physical parameters correspond to those 
in Sakai et a/. [11]. Fig. [1] shows the time history of magnetic field energy, ion and 
electron total energies, respectively. As discussed in Sakai et al. [H] , the electron energy 
is enhanced by the excitation of MTSI (the arrow in Fig. [1]). Fig. [2](a) clearly shows 
the obliquely propagating waves with parallel electric field, corresponding to Fig. 3(b) 
in Sakai et al. [11] . We have also confirmed that the propagation angle of the waves is 
about 60°, which agrees with the simulation results in Sakai et a/. [H] and the linear 
analysis in Wu et a/. [21] (We note that while Wu et a/. [21] discuss the case c^e >> l^el. 
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Figure 2. Snapshots of longitudinal electric field Cx at {a)t — 32rjj ^ and (b)i = 64rij ^ 
in Run f . 



the growth rate of the MTSI is not sensitive to co'^/|n^|P^).The initial perpendicular 
velocity (Vq) is obtained from eq.(jl]) 

(l + r)a;o 



V, 



[Uio - MeOj 



n, 



^ (5) 

where v^o = uo/ko, r = nii/me. Vq in Run 1 is Vq = 0.355Ca. We remark that the 
kinetic effects of electrons are important for the MTSI observed here, since Vth,e > Vn|21j. 
After these obliquely propagating waves excited by the MTSI are damped by 
electron Landau damping, low-frequency parallel propagating waves become dominant 
(Fig. [21(b)). Dissipation of these parallel propagating ion acoustic waves heat the ions 
as discussed in Sakai et a/. [14]. Although Sakai et a/.[T3j argued that such parallel 
propagating longitudinal waves are excited by transverse modulation of parent Alfven 
wave, we consider these waves are more likely to be generated by the usual, parallel decay 
instability. Fig. [3] shows the power spectrum of the complex magnetic field b (= by + ibz) 
in the k^ — ky wave number space at t = 96f2~^. The anti-parallel propagating daughter 
waves {ujo/ko < 0) are dominantly excited around {kx,ky) = (0.73,0). Presence of 
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Figure 3. The power spectrum of the complex magnetic field b (= by + ibz) in kx — ky 
wave number space at i = 96ri^ in Run 1. The vertical and horizontal axes {ky and 
kx) are normalized by the ion skin depth. 



such anti-parallel propagating waves provides a strong evidence of the occurrence of 
the decay instability, while the daughter Alfven waves excited by transverse modulation 
have finite propagation angles. In addition, these daughter waves are consistent with 
the linear analysis of Hall-MHD system for parallel propagating waves p^ ITS] . Thus, 
we conclude that under this set of parameters, the one-dimensional decay instability 
dominantly dissipates the parent Alfven wave, as the Hall-MHD theory predicts [TSj. 
While numerical results of the present study and Sakai et al. [S] are quite similar, their 
interpretation on the parametric instability of parent Alfven wave seems misleading. 

Next, we discuss the RH- mode case (Run 2). Fig. H] shows the time history of 
energies in the same format as Fig. [TJ In contrast to Fig. [1] (Run 1), electron heating 
due to the excitation of MTSI is not evident. From eq. (EI), Vq in Run 2 is 0.223Ca, 
which is larger than Vth,i- We here find that even when Vq > Vth,ii which is a rough 
criterion of the MTSI[23], the cross field current is not large enough to excite the MTSI 
under the particular set of parameters used. On the other hand, the decay instability 
takes place in a way similar to Run 1 (not shown). 

In Run 3, we discuss the case in which the Walen relation is not satisfied for the 
initial Alfven wave. These initial conditions are intended to simulate the plasma and 
the magnetic field fluctuations in the solar wind and solar surface, where both RH- and 
LH- polarized Alfven waves are considered to be mixed|2]. In such a case, amphtude 
of transverse magnetic field and transverse bulk velocities vary in time|20j. Since such 
linear oscillations modify the magnitude of the cross field current, they let the system 
stay in the stable and unstable regimes in turn, so that the MTSI growth is effectively 
reduced. 

Finally, we discuss the results of Run 4, in which parameters are the same as those 
in Run 1 except for vth,i = 0.1c. In contrast to Run 1, MTSI is suppressed by the 
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Figure 4. Time liistories of magnetic field energy (thick line) , ion (dashed line) and 
electron (thin line) total energies in Run 2. 
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Figure 5. The power spectrum of the complex magnetic field b (= by +ibz) in k^ — ky 
wave number space at i = 160fl^ in Run 3. The vertical and horizontal axes {ky and 
kx) are normalized by the ion skin depth. 



ion kinetic effects (Fig. [6]). Further, as Fig. [7] shows, the decay instabihty is not the 
dominant instabihty in this case, because of the occurrence of "kinetic" modulational 
instabihty [20]. In Fig. [71 both the dominant daughter waves {k^ ~ —0.25,-1.25) and 
the higher harmonic waves are observed. 



4. Summary 

In this paper, we discussed the parametric instabilities of Alfven waves using a 2-D PIC 
simulation code. First, we re-ran the simulation of Sakai et a/.[Ilj using essentially the 
same set of parameters. While our numerical results are basically the same, we conclude 
that the generation of longitudinal waves are due to the usual, parallel decay instability, 
rather than the transverse modulational instability proposed by Sakai et a/.[ll], since the 
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Figure 6. Time liistories of magnetic field energy (thick; line), ion (dashed line) and 
electron (thin line) total energies in Run 4. 
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Figure 7. The power spectrum of the complex magnetic field b (= by +ibz) in fca^ — ky 
wave number space at i = 128r2~ in Run 4. The vertical and horizontal axes (ky and 
kx) are normalized by the ion skin depth. 



one-dimensional decay instability is dominant in a low beta plasma as the Hall-MHD 
theory predicted [TOl [T3] . Further, we performed several other simulation runs using 
different parent Alfven wave polarizations and different ion thermal velocities. The 
condition to excite the MTSI is easily violated by such modifications of the parameters. 
Thus, the electron heating by the MTSI of the Alfven waves seems unimportant in 
the upper chromoshpere, since the Alfven waves excited by the convection in the solar 
surface are likely to be the superposition of Alfven waves with different polarizations 
and propagating directions, but not the monochromatic, LH- polarized wave. Further, a 
similar process on obliquely propagating waves have also been suggested by Markovskii 
and HoUweg [2^ according to the proton heating. The more detailed analysis on this 
issue is desired. 

While the dominant parametric instabilities are quasi-one dimensional in a low 
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beta case, the transverse modulational instability can be more important in different 
circumstances. For instance, the Alfven wave excited by the ion beam instabilities 
can be dominantly dissipated through the transverse modulational instability, which 
is nonlinearly driven by the obliquely propagating waves excited by the same beam 
instabilities [25]. Further, even if the parametric instabilities are one dimensional, two 
or three dimensional density structures are excited as a consequence of its nonlinear 
evolution im [26] . When the inhomoginity of the background plasmas exist, heating of 
electrons by phase mixed Alfven waves take place [27]. We remark that even in the 
uniform background plasmas, the electron heating regardless of the MTSI is actually 
observed in PIC simulation {e.g., Sakai et a/. [Hj and Fig. Hj). The electron heating 
is possibly caused by numerical noise originated from the super-particles in the PIC 
simulation. To avoid this and confirm the recent results, the electromagnetic Vlasov 
simulations are in plan. 
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